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ABSTRACT 

We report a ~ 3a detection of a 21.82 A absorption feature in the direction 
of the BL Lac MRK 421 which is interpreted as O VII Ka at z ~ 0.01. This 
corresponds to the redshift of a H absorber in a cosmic void detected by HST 
along the line of sight of the same object. The 21.82 A line proves the existence of 
warm/hot intergalactic medium (WHIM) outside the Local Group, in agreement 
with current models of structure formation in the Universe. The WHIM at 
z ~ 0.01 has temperature of z (2.5 — 4) x 10 6 K for gas densities in the range 
10~ 6 - 1 atom cm . We also detect O VII, O vm and Ne IX absorption lines 
with zero velocity, possibly connected with WHIM in the Local Group. 

Subject headings: BL Lacertae objects: individual (MRK 421) — intergalactic 
medium — large-scale structure of universe — quasars: absorption lines 

1. Introduction 

In the local Universe ~ 50 — 60% of the baryons visible at z > 2 and predicted by 
the Standard Big-Bang nucleosynthesis are undetected (e.g. Fukugita, Hogan, & Peebles 
1998). Numerical simulations predict that such baryons, in the form of Lyct clouds (i.e. 
gas at T < 10 4 K) and visible in the UV spectra of bright AGNs, have been shock-heated 
during the formation of the structures in the Universe to temperatures of ~ 10 5 — 10 7 K (i.e. 
the warm hot intergalactic medium, WHIM). The most efficient way to detect such gas is 
through resonant absorption lines from highly ionized metals (e.g. 0vi,0vil,0vill,Nelx) in 
the far UV and soft X-ray spectra of background sources (e.g. Hellsten, Gnedin, & Miralda- 
Escude 1998; Fang & Canizares 2000). Current far UV observations proved the existence 
of the low temperature tail (T ~ (1 — 5) x 10 5 K) of the WHIM through the detection of 
Ovi up to z ~ 0.2 (e.g. Sembach et al. 2000; Tripp et al. 2001). However the bulk of the 
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WHIM baryons (~ 70%, e.g. Fukugita, Hogan, & Peebles 1998) are expected to lie at higher 
temperatures (T ~ 10 6 — 10 7 K) and their features should be detectable in the soft X-ray 
regime. Chandra and XMM-Newton high resolution spectroscopy (R ~ 100 — 500) brought 
to the detection of a "zero redshift" Ovn absorption feature at 21.6 A in the spectra of four 
bright AGNs (PKS 2155-304, Nicastro et al. 2002; Fang et al. 2002; Cagnoni et al. 2003; 
3C 273, Fang, Sembach & Canizares 2003; H1821+643, Mathur, Weimberg & Chen 2002 and 
MRK 421, Nicastro et al. 2001), interpreted as the signature of WHIM present within the 
Local Group. The only published evidence of WHIM outside the Local Group is at z ~ 0.05 
(Fang et al. 2002), however this feature is not confirmed in other Chandra and XMM-Newton 
observations (e.g. Nicastro et al. 2002; Cagnoni et al. 2003). In this paper I present strong 
evidence for an Ovn Ka absorption line at z ~ 0.01 in the spectrum of the BL Lac object 
MRK 421. 

The paper is organized as follows: § 2 reports the XMM-Newton observations of MRK 421 
and describes the data reduction; § 3 contains the spectral fits and the discussion on the 
observed WHIM absorption features. § 4 is a conclusive section containing a summary of 
the results. 



2. Observations and data reduction 

MRK 421 is a calibration target for XMM-Newton and has been observed several times 
from the launch (Dec. 10, 1999) up to the time of writing. Table 1 summarizes all the 
public observations as of December 3, 2002. MRK 421 X-ray spectrum is thought to be 
relativistically beamed synchrotron emission from energetic electrons and its intrinsic lack of 
features is the ideal laboratory to search for faint WHIM absorption lines. In this paper I will 
concentrate on the high resolution (AE/E from 100 to 500, FWHM, or 100 to 800, HEW, in 
the energy range 0.33-2.5 keV - 5-38 A) Reflection Grating Spectrometers (RGS1 and RGS2) 
data collected by the two XMM-Newton X-ray telescopes 1 . XMM-Newton RGS effective 
areas are complex in shape, and contain tens of narrow dips due to bad or hot columns or 
pixels in the CCD detectors, which require extremely accurate calibration measurements for 
a proper modeling (e.g. Cagnoni et al. 2003). Current RGS calibration uncertainties are as 
accurate as ~ 5 — 10% between 7 and 36 A and, as a consequence, false absorption/emission 
features with such relative intensities, are expected in the RGS spectra in physical units, 
in proximity of the known instrumental features. The strongest resonant absorption lines 
from neutral and/or highly ionized O and Ne, fall in wavelength ranges (i.e. 13-14 A, 18-20 



1 I refer the reader to Scmbay et al. (2002) for a detailed description of the observations and for and 
accurate spectral and timing analysis of EPIC data collected in 2000 and 2001. 
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A, 20-24 A) in which RGS-2 spectra either do not exist (20-24 A, due to the failure of a 
CCD chip) or contain strong line-like shaped instrumental features (see Cagnoni et al. 2003). 
Therefore I rely on RGS-1 spectra, which, instead are relatively instrumental- feature- free in 
these wavelength ranges, and use RGS-2 to double check the reality of a line, when possible, 
and to cover the 10.5-14.2 A region, where RGS-1 has a failed CCD chip. I also restrict our 
analysis to the first order spectra only. 

I reprocessed the data using XMM-Newton Science Analysis System (SAS) version 5.3.0 
and the latest calibration files as of December 3, 2002. Since the wavelength calibration of 
XMM grating spectra strongly depends on the position of the 0th order, I used the VLBI 
position as centroid of the 0th order source (Ma et al. 1998). Extraction regions, for source 
and background, were chosen to be, respectively, within the 95% and outside the 98% of 
the PSF. To exclude high particle background periods caused by solar activity, I extracted 
the background lightcurves from CCD-9 and excluded all the time intervals for which the 
background count rate was higher than 1.5 x 10~ 3 count s -1 . The net exposures for each 
observation are reported in Table 1. 

In order to improve the signal to noise ratio (SNR) I combined all the RGS-1 order 1 
spectra obtained with a pointing offset < 15" 2 (see Table 1) and computed the combined 
RGS-1 response matrixes (a convolution of the detector response matrix and effective area, 
Fig. 1) using the PINTofALE interactive data language software suite (Kashyap & Drake 
2000). I did the same for RGS-2. The total net exposures are 125 ks and 122 ks for RGS-1 
and RGS-2 respectively. 

3. The WHIM detection 

I fit the combined fluxed RGS-1 spectrum in the wavelength range 14.2-38.1 A using 
version 2.2 of SHERPA (Siemiginowska et al., in prep.) modeling and fitting tool from the 
CXC analysis package CIAO 2.2 (Elvis et al., in prep.). I used an absorbed power law 
model with absorption characterized by a column density N(X) and an absorption cross 
section <7x for each element. I included in the model H, He I and He II (Rumph, Bowyer, 
& Vennes 1994) and heavier elements (Morrison & McCammon 1983). I fixed the Galactic 



2 The inclusion of the off- axis data in the combined spectrum increases the number of photons at 20 A 
from ~ 4000 to <~ 7000 per 0.06 A resolution element, but degrades the quality of the spectrum. May 2002 
observations, for example, are highly affected by hot columns and pixels which result in spurious features in 
and around the expected WHIM lines. To reduce the impact of hot columns and pixels, the RGS units have 
recently been cooled (November 2002). 
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hydrogen column density at N H = 1.61 x 10 20 cm 2 (Lockman & Savage 1995) and the ratios 
N Hel /N m =0.1 and N HeU /N m =0.01. The best fit RGS-1 photon slope is T A = 0.095 ±0.004 
(corresponding to an energy photon slope of 2.095 ± 0.004) and the normalization is 0.0136 
photons cm -2 s _1 A -1 at 20.5 A. In order to properly model the continuum around the 
expected WHIM features, I performed local absorbed power-law fits in the 18.0-20.6 A and 
m the 20.9-22.4 A regions. I kept the absorbing column fixed to the Galactic value in 
MRK 421 direction and modeled the absorption lines with Gaussian. The lines detected at 
> 2a and the corresponding best fit parameters are listed in Table 2. 

It appears to us that the only possible interpretation of the 21.82 A line (2.8a detection; 
Fig. 1) is redshifted O VII Ka. This places the gas producing it at z ~ 0.01 and makes this 
a solid X-ray detection of WHIM outside our local group of Galaxies. The corresponding 
z ~ 0.01 O VII K(3 is expected to be too faint (EW~ 0.8 mA at ~ 18.84 A) to be detected in 
this spectrum. If any O Vlll were present in the z ~ 0.01 WHIM, the strongest expected line 
would be Lya at 19.18 A, not detected in our spectrum. The 3a upper limit on the line EW, 
computed fixing the line FWHM to the 21.82 A line FWHM value, is 2.47 mA. Under the 
assumption of unsaturated line the EW ratios between different species depend on the gas 
temperature and density. Using the ratio between the upper limit on the EW of O vm at 
z ~ 0.01 and the 21.82 A line EW, I obtain an upper limit on the z ~ 0.01 gas temperature 
of ~ 4 x 10 6 K, for a gas density of 1 atom cm -3 and of ~ 2.5 x 10 6 K, for a gas density 
of 10~ 6 atom cm -3 (see Fig. 5 in Nicastro et al. 2002). These values are consistent with 
the temperature range predicted for the WHIM emitting in the soft X-ray band. I derive 
from the 21.82 A line a velocity of cz ~ 3062 km s _1 , consistent, within the errors, with the 
redshift of a strong H Lya line detected by HST on MRK 421 spectrum {cz = 3035 ± 6 km 
s _1 ) and coming from absorbing material in a cosmic void (Shull, Stocke, & Penton 1996; 
Penton, Stocke, & Shull 2000). The location of WHIM in a cosmic void is consistent with 
the picture that such gas is connecting the overdense regions of the sky that collapsed into 
clusters and groups during the structures formation. 

For the other lines listed in Table 2, the strongest is the O VII Ka absorption at 21.6 A 
(Fig. 1); the corresponding O VII K(3 is visible at 18.67 A (Fig. 2). The 21.6 A feature has 
already been seen in MRK 421 itself (Nicastro et al. 2001) and in other 4 AGNs (Nicastro 
et al. 2002; Fang et al. 2002; Cagnoni et al. 2003; Mathur, Weimberg & Chen 2002; Fang, 
Sembach & Canizares 2003) and it is usually attributed to a WHIM within our local group 
of galaxies (e.g. Nicastro et al. 2002) or to radiatively cooling gas inside our Galaxy (e.g. 
Heckman et al. 2002). I also report the detection of zero redshift O vm Lya at 18.97 A 
(Fig. 2) and of Neix Ka at 13.4 A (Fig. 3). Note however that these lines are affected by 
small features in the effective area. 

The EW of the zero redshift features (O vn,0 Vlll and Neix) are compatible with those 
measured with Chandra for PKS2155-304 (Nicastro et al. 2002; Fang et al. 2002) and for 
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MRK 421 (Nicastro et al. 2001). In particular, the EW derived fixing the lines FWHM to 
the value of the 21.6 A line, (i.e. 0.33, 0.27 and 0.20 eV respectively) are consistent with 
the lines being unsaturated and with column densities of ~ 10 16 atoms cm -2 per ion species 
(see Fig. 4 in Nicastro et al. 2002). For a gas density of ~ 10~ 6 atoms cm -3 , Neix/Ovil and 

vm/Ovil ratios indicate gas temperatures of ~ (0.4 - 2) x 10 6 K, and ~ (1 - 4)10 6 K, 
respectively (Fig. 5 of Nicastro et al. 2002). These ranges move to ~ (1.5 — 2.5) x 10 6 K and 
~ (3 — 4) x 10 6 K in case of Galactic gas density. Even if these estimates have to be regarded 
with caution because of the possible modifications of the effective area features, they are all 
consistent with the range of temperatures predicted for the WHIM. 

4. Conclusion 

I report strong evidence of WHIM outside our local group of galaxies through the de- 
tection of a 21.82 A O VII Ka absorption line in the XMM-Newton RGS-1 spectrum of 
MRK 421. The line redshift is ~ 0.01, in agreement with that derived for the strong H Lya 
absorption line detected by HST and related to an absorber in a cosmic void. The three 
closest galaxies are at 2.14, 3.98 and 4.12 Mpc (Penton, Stocke, & Shull 2000). This 
detection is a firm proof the existence of a WHIM filament at z ~ 0.01 possibly connecting 
the nearby galaxies with other denser regions of the sky, as expected by the models of the 
cosmic structures formation and evolution. I derive an upper limit on the gas temperature 
of ~ 2.5 — 4 x 10 6 K for a gas density of 10~ 6 — 1 atom cm -3 . The only known similar 
evidences are a 4.5a detection of O vm at z ~ 0.05 reported by Fang et al. (2002) in the 
Chandra spectrum of PKS 2155-304 and the ~ 2a detections of two O VII lines at z ~ 0.2 
in the Chandra spectrum of H1821+643 (Mathur, Weimberg & Chen 2002). However Fang 
et al. (2002) feature is not confirmed in other Chandra and XMM-Newton observations (e.g. 
Nicastro et al. 2002; Cagnoni et al. 2003). 

1 also detect the zero redshift absorption features of O VII, Ovil and Neix with positions and 
EWs consistent with previous detections in the direction of other bright AGNs (e.g. Nicastro 
et al. 2002; Cagnoni et al. 2003). 



I acknowledge a C.N. A. A. fellowship. I thank Aldo Treves and Francesco Haardt for 
useful scientific discussions and a careful reading of the manuscript. 
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Fig. 1. — Top panel: MRK 421 20.9-22.6 A combined RGS-1 spectrum and best fit absorbed 
power law. Middle panel: residuals to the fit and, bottom panel, RGS-1 combined effective 
area. 
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Fig. 2. — Same as Fig. 1 for the 18.0-19.45 A region. 
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Fig. 3. — Same as Fig. 1 for the RGS-2 vawelength range 12.0-14.2 A. 
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Table I. All the public XMM-Newton RGS observations of MRK 421. 



Obs. Id. 


Obs. Start time a 


Total Exposure 13 


Net Exposure 


Offset 


Used observartion? d 


0099280201 


00-11-01 


40 


34 


6.3 


YES 


0099280301 


00-11-13 


50 


40 


5.87 


YES 


0099280401 


00-11-14 


43 


27 


116.69 


NO 


0099280501 


00-11-13 


21 


15 


7.63 


YES 


0099280601 


00-11-15 


20 


12 


117.4 


NO 


0136540101 


01-05-08 


39 


36 


12.95 


YES 


0136540201 


01-05-08 


10 





13.16 


NO 


0153950601 


02-05-04 


40 


39 


123.98 


NO 


0153950701 


02-05-05 


20 


18 


169.43 


NO 


0153950801 


02-05-05 


22 


21 


177.49 


NO 



a YY-MM-DD 

b Total RGS-1 (-RGS-2) on time in ks. 

c For RGS-1 in ks after the high particle background times rejection. RGS-2 net exposure times 
are similar. 

d For details see footnote 2 in the text. 
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Table 2: Best-fitting RGS-1 absorption line parameters and la errors. 



Line ID 


X 




FWHM 


EW a 




a 




(A) 


(km s _1 ) 


(io- 2 A) 


(mA) 






Neix b 


13.441 ±0.029 


-163±g» 


12.061^3 


3.52+1;™ (2.89li;») 


2.5 


Ovn K(3 


18.6591H 




i i 09+2.8O 
ii.zz_ 10 _o7 


5.47^ 6 f (4.88; 


fl.29\ 
-1.37) 


3.8 


OviII Lya h 


18.961 ±0.008 


-98+ 41 


o 97+IO.OI 
-17.34 


5.50+f^ (7.98^ 


-1.31\ 
-1.33^ 


4.8 


Ovil Ka 


21.597 ±0.008 


-75 ± 80 


8.77 + }-| 


12-67±i:f2 




7.6 


Ovil Ka 


91 oon+0.027 
Zl.OZZ_ 024 


3062+j| 


2.77l?;l 


3.681*:™ 




2.8 



a The values in parenthesis are obtained fixing the FWHM of the line to the value of the zero rcdshift Ovn 
Ka 

& These lines could possibly be enhanced or modified by the presence of small effective area features (see the 
bottom panels in Fig. 2 and Fig. 3) 



